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j^ ■ Abstract 

We consider a non-linear, one-dimensional wave equation system with finite-dimensional 

stochastic driving terms and with weak dissipation. A stationary process that solves 

the system is used to model steady-state non-equilibrium heat flow through a non- linear 

medium. We show existence and uniqueness of invariant measures for the system mod- 

i-^ ' ified with ultraviolet cutoffs, and we obtain estimates for the field covariances with 

respect to these measures, estimates that are uniform in the cutoffs. Finally, we discuss 

,^ ■ the limit of these measures as the ultraviolet cutoffs are removed. 
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^ ■ 1 Introduction 

O 

In this article, we consider a dynamical system of equations for a scalar field in one space 

X/^ • dimension with ultraviolet cutoff, where the equations have a non-linear term, dissipative 

^^ I terms, and stochastic driving terms. We will be concerned with the problems of constructing 

^-^ ' measures on the space of field configurations invariant under the time evolution of the 

system, and of estimating the covariance of the field with respect to these measures. 

This system of equations provides a model for heat flow; the system is a variant of 
rS I models introduced and examined in detail by Eckmann, Pillet, and Rey-Bellet ([10l[Tl], see 

C^ ■ also [Ml ESI [Ml [3 [22] ) . These authors considered a finite system of non-linear oscillators 

coupled to two (or more) free fields that are governed by linear wave equations. The free 
fields, which model heat reservoirs, are given Gaussian-distributed random initial conditions 
(Gibbs states), but in general at different temperatures Tj for each of the free fields. The 
equations of motion for the free fields are readily integrated, leaving a system of stochastic 
equations for the oscillators. For suitable couplings, these authors found that the equations 
were Markovian stochastic differential equations. The equations given below, then, are sim- 
ply the analogue of these stochastic differential equations, but with the oscillators replaced 
by a non- linear scalar field in one space dimension with bounded support. We examine 
these equations with an ultraviolet cutoff, i.e. where Fourier modes of the field are set to 
zero for mode numbers {n}, n > M, and we consider the limiting situation where the cutoff 
M is reixLoved as M — )• oo. 
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The equations of motion for the system are given by 

dt4>{x,t) = 7r(x,t) 

dt7r{x,t) = {dl — l)(f>{x,t) — g{(f){x,t)) — r{t)a{x) 

dr{t) = -{r{t)-{a,TT{t)))dt + Vfduj{t). (1.1) 

Here (^, 7r)(x, t) is a pair of scalar fields; we assume, for convenience, that they satisfy peri- 
odic boundary conditions for x G [0, 27r]. The fixed vector- valued function (or distribution) 
a = (ai(x), 02(2;)) has two components, each in a Sobolev space Hg, that will be further 
specified below. The (dependent) random vector r{t) = {ri{t),r2{t)) takes values in M^; its 
components are artifacts of the reservoirs. The Oj couple the oscillators to the reservoirs, and 
here r{t)a = Yli=i i"i{'t)'^i{x)- The quantity (a^, vr) is simply the i^([0, 27r], drc)-inner prod- 
uct of ai{x) and 'ir{x,t). The driving term uj{t) = {uJi,uJ2){t) is standard two-dimensional 
Brownian motion and the Tj, i = 1,2, are the reservoir temperatures. The non-linearity g 
is assumed to be an odd, bounded, Lipschitz-continuous function. 

The overarching goal of the study of these equations is to construct a stationary (non- 
equilibrium) process $(t) = {cj}{t),Tr{t),r{t)) governed by the above equations and having 
physically reasonable sample path properties. In particular, the field component (p = (j){x, t) 
should be at least square integrable in x a.s., or actually continuous in x a.s. as in the 
equilibrium case Ti = T2 and in the non-equilibrium but linear g = case, all of which we 
review below. 

In this article, we establish the existence and uniqueness of stationary states for these 
equations with ultraviolet cutoff M, i.e. such that Fourier modes of the field (p and its 
momentum vr are set to zero for mode numbers |n| > M. We also obtain estimates on the 
field covariance uniform in the cutoff with respect to the stationary cutoff measures. For 
certain choices of the coupling functions a, we show that (/> is indeed L^-integrable a.s., 
and that ||(/'||2 has variance uniformly bounded in M. We believe that these estimates are 
an important step in establishing bounds on the mode variances uniform in both the mode 
number n and the cutoff M, as one might expect on physical grounds, an issue that we 
will address in subsequent work. We conclude with a discussion of the limit of the cutoff 
measures, M — >• 00. 

It is known that in equilibrium, temperatures Ti = T2, and even for a non- focusing 
unbounded nonlinearity g{x) oc x^, there is an invariant measure for the system (essentially a 
Gibbs state) which is supported on field configurations {i;^} that are almost surely continuous 
in X (they are Brownian motion-like, in fact. Holder continuous with any index < 1/2). The 
measure is absolutely continuous with respect to the g = Gaussian measure (see p7]). 
For the non-equilibrium case Ti 7^ T2 but with g = 0, there is steady-state energy flow, 
i.e., there is a random variable measuring heat flow through any fixed point x G [0, 2tt], and 
with respect to this measure, the expectation of this heat flow is non-zero. Moreover, the 
measure has support properties similar to those of the equilibrium case (see [30]). But for 
the non-linear 5 7^ non-equilibrium case, the regularity of the field is at present not known. 
If the above equations are to serve as a physical model for a non-linear vibrating string in a 
non-equilibrium stationary state, we expect the measure to have support properties on field 
configurations similar to those in the equilibrium case. In particular, we expect the field 
configurations of to be a.s. continuous, not shattered. It is really this latter problem- 



establishing the regularity of the field with respect to the invariant measure- that is our 
primary concern in this article and its sequel. The first step in this program is to obtain 
estimates on the field covariance uniform in the ultraviolet cutoff. 

We thus consider ultraviolet cutoff versions of the system (jl.ip where we retain the 
Fourier modes {(f){n)} for the field (p and {7r(n)} for momentum vr for |n| < M with M a 
positive integer. Let P<m be projection onto the Fourier modes {n}, with \n\ < M. Then 
setting $M = {(pM^T^M^fM)'^ as the cutoff field having non-zero modes only for |?7,| < M, 
otMi = PmoHj we have that our equations can be written in matrix form as 



d^M = (Am^m - Gm($m)) dt + VTdw, 



where 
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with {aM,i the linear functional of integration against aM,i-, and with 

/ \ / ( 
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We make the following assumptions about the coupling functions and the nonlinearity g for 
the system of equations (jl.ip and ()1.2p : 

Assumption 1.1. The coupling functions Oi, i = 1,2 are real, with all of their Fourier 
coefficients {ai{n)} non-zero. There exists a positive constant ci (independent of n) such 
that 

ci(|ai(n)|2 + |a2(ra)P) < \ai{n)'^ + Mnfl, (1-5) 

and there exist both a positive constant C2 and 0, —1/2 < 9 < 1/2, such that, for all n, 

C2{n' + lf'<\a,{n)\<-{n' + lf/\ (1.6) 

C2 

The non-linearity g is assumed to be bounded and uniformly Lipshitz. 

The assumptions on the coupling functions are made in order to assure that the per- 
turbed eigenvalues of A and Am are non-degenerate and have negative real parts. In our 



model, the dissipation of large n modes is very weak, with rate 



-n 



29-2 



It is thus impor- 



tant that the coupling functions themselves have Fourier coefficients large enough so that 
the dissipation is not overwhelmed by the noise driving these modes, and at the same time 
not so large that the equations of motion make no sense as the ultraviolet cutoff is removed. 
The assumptions on the a's above appear to capture the right growth/decay rates in the 
a(n)'s. 



The problem of constructing an invariant measure for this gently stochastic wave equa- 
tion may be contrasted with the problem of constructing stationary states for stochas- 
tic parabolic and even 2-dimensional stochastic Navier-Stokes equations with viscosity, 
where the dissipation of large n Fourier modes of the field is strong (for Navier-Stokes, 
see (U [20l \12\ HH IE])- For a heat equation the n*'^ mode decays at a rate ~ — n^, com- 
pared with our rate ~ — n^^~^ mentioned above. In [1], an invariant measure is obtained 
for a non-linear wave equation with cylindrical Brownian motion driving terms but with 
strong dissipation for all Fourier modes. In |27j , an invariant measure is constructed for a 
one-dimensional stochastic non-linear Klein-Gordon field but at thermal equilibrium, with 
weak dissipation at high Fourier modes; see also [3^ where an invariant measure for the 
non-equilibrium but linear case is given as an invariant measure for an Ornstein-Uhlenbeck 
process. See also [18], for an invariant measure for a harmonic crystal in non-equilibrium. 
For invariant measures for wave equations in equilibrium but with no stochastic driving 
terms, i.e., equilibrium statistical mechanics for wave equations, see for example [21^ 131]. 
and for the non-linear Schrodinger equation, see pt [5t[T9]. 

In the following section, we show existence and uniqueness of invariant measures for 
the ultraviolet cutoff systems, and in section 3 we obtain bounds on the field covariance 
modes with respect to these invariant measures. In section 4, we give a general discussion 
of tightness for these cutoff measures. Much of our analysis will use perturbation theory 
for linear operators, with estimates on the tails of the relevant perturbation series. The 
appendix summarizes these perturbation calculations. 

2 Existence and Uniqueness of the Stationary Measures for 
Systems with Ultraviolet Cutoff 

In the following, let Ti be the Hilbert space of complex-valued functions {<!> = ((/), vr, r) } 
equipped with inner product 

{^i,^2)h= id^(l)A(l)*2 + (j)icl)*2+7riiT*2)dx + ri-r*2, (2.1) 

Jo 

and let || • ||-^ denote the corresponding norm for Ti. 

The solution $M(t) = ^M{t; *) to Equation (fL2]) with initial data ^m{0) = * satisfies 
the Duhamel integral equation 

$A/(t) = - [ e^'~'^^^'GM{'^M{s))ds+ [ e^'-'^^^'VTdu{s) 
Jo Jo 

+ e*^«*, (2.2) 

(see [29], chapter 5). The field <&Af(0 is a vector in the (4M + 4)-dimensional subspace 
T-Lm of Ti spanned by (2M + 1) Fourier modes each for (/) and vr, and by the two r modes. 
Using <I>Af(t), one defines a probability semigroup P^j acting on bounded Borel functions 
on m4A^+4 via 

Pijf{^) = E[f{1>Mit;^))] (2.3) 



with E[-] the Brownian motion expectation. The operator P^j is a Feller semigroup, i.e., 
P^^jf{^) is continuous in ^ if / is continuous which follows from the integral equation above 
for ^M{t) and the assumption that g is Lipschitz. 

Proposition 2.1. There exists an invariant measure for the semigroup P\,j. 

Proof: To show the existence of an invariant measure it is sufficient to prove that the family 
of measures {^^m^t} defined by 



fdfiM,T^^I Plifmdt 



(2.4) 



T>0 



is tight on the space T-Lm [Sj, for /'s bounded and continuous. Since our system is finite 
dimensional, one need only show that for each ^ G T-Lm^ and e > 0, there exist an i? > 
such that for time T sufficiently large. 



1 

T 



T 



P{\\^M{t-,mH>R)dt<e. 



(2.5) 



By the Markov inequality, we have 



P (||$Af (t; mu >R)<^E[ \\^M{t- mn] , (2.6) 

so it suffices to show that E [||$jv/(i) ^)II'h] is bounded uniformly in t. But from ()2.2p . 



£;[||$M(i;*)|y < E 

+ E 



I e(*-^)^^^(GM($M(5)))rf 
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I e(*-")^«Vf(ia;(s 
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(2.7) 



Now the unperturbed matrix operator Bm obtained from Am by setting the coupling 
a's to zero has spectrum consisting of a doubly degenerate eigenvalue A_i(0) = —1, simple 
eigenvalues Aq (0) = ±z, and doubly degenerate eigenvalues A^(0) = ±i\/n'^ + 1, 1 < n < 
M. As seen in perturbation theory, the matrix Am has eigenvalues near those of 'Bm, but 
with small negative real parts (see the Appendix, Lemma (lA.ip . Eq. (IA.7p ). 

Given these eigenvalue estimates, it is then easy to see that the first and third terms 
on the right side of (12. 7p are bounded (Gm is bounded) uniformly in t. The second term 
involving the Ito integral is estimated via Ito calculus. 



E 



z^^-''^^^'VTduj(s] 



n 



< E 



,{t-s)Aj, 



Tdoj{s] 



n 



Trn I e' 

10 



s Am rr^s aI 



'Te'^^Mds, 
which is uniformly bounded in t, again by the eigenvalue estimates for Am- 



{2.i 



Remark: The above argument can also be used to show boundedness of higher moments of 
the field $m with respect to the invariant measure. 

We next proceed to showing uniqueness of the invariant measure fiM satisfying (imPm — 
fiM for the semigroup P|,j. 

Proposition 2.2. Under the assumptions il-l\) on the coupling functions a and nonlinear- 
ity g, the invariant measure ^m for P\^ is unique. 

Proof: In order to show uniqueness of the invariant measure, it suffices to show that the 
system in which the random field ^Kiif) is replaced by a deterministic function Xm(0 G T-Lm 
satisfying 

dXM{t) = {AMXM{t) - GM{XM{t)))dt + VTu{t)dt (2.9) 

is controllable. This means that, given Xq, Xi, and a time ti, one can find a smooth 
control u{t) = (ui{t),U2{t)) such that Xjv/(i) is a solution to this differential equation, 
and Xjvf(O) = Xq, XMih) = Xi (see [6j, page 144). But controllability of this system is 
equivalent to controllability of the simpler system 

dXM{t) = (AMXMit) - GMiXMit)))dt + VTuit)dt, (2.10) 



where the drM(t)- equations are replaced by drM = vTu{t)dt, i.e., the new and the old 
controls are related by 



Tn' = - {rM{t) - /3(aAf , vrM(t))) + VTu(t), (2.11) 

and the modified matrix Am is obtained from Am by replacing its bottom two rows by 
zero rows. The solution of equation (|2.1U|) satisfies 

XMit) = - [e^'-'^^^^GM{XM{s))ds + [ e^'-'^^^' VTu{s)ds + c'^^Xq. (2.12) 

Jo Jo 

We will show the controllability by first replacing Xjvf(s) in the first integral on the 
right side of this last equation ()2.12p by an arbitrary continuous function Zm^s), thereby 
obtaining an expression for u{t) explicitly in terms of ZAf(s), and then by using a fixed 
point theorem. We thus consider 

XM{t) = - [e^'-'^^^'GM{ZM{s))ds + [e^'-''^^^'VTu{s)ds + e*^«Xo, (2.13) 

Jo Jo 

for which we seek a control u{t) such that for fixed Zm{s)^ XM{ti) = Xi. 

Let 

/■* - -t 

Dt = / e-'^'''Te-'^Mds. (2.14) 

Jo 



Lemma 2.3. The matrix Dj is a non-negative self-adjoint invertible (4M + 4) x (4M + 4) 
matrix for t > 0, provided the Fourier coefficients of ai and 02 satisfy the assumption 
Eq.mMl- 



Proof: Clearly, the integrand of the integral in (j2.14p is non-negative and continuous in 
s. If Dt were not invertible, then for some X ^ , {X, D^X)-^ = 0, and in partic- 
ular VTc'^^-^mX = for all < s < t. Since a/Tc^^^mX is smooth in s, we have 
^ (VTc-'^mx) ^ = for ah n, so ^/TAj^"X = 0, n = 0, 1, 2, • • • . However we know 

t ^ 

{{AMrVTY,X)n=0 (2.15) 



that for any Y £ TIm, {Y, VT{A\j"-X)n = 0, which imphes that 



for all Y, and thus the range of {(AAf)"vT}^Q is not the whole space Hm- But this 
contradicts the following claim. 



Claim 2.4. Vectors consisting of the columns o/{(A^/)'^vT}5^o n = 0,1,2, .. . span Hm- 

Remark: The generator of the Ornstein-Uhlenbeck semigroup defined by the linear equation 

d<^M{t) = AM^M{t)dt + y/T(hj is 

d^ d^ f 6 



' drj dr^ ' J "' ScpM 

dx {{dl - l)(t>M + ^rau) j^ ■ (2.16) 



If we denote 



Xi = -— , X2 = -— , and 
or I or 2 

Y = j dxLM^ + {{dl-l)<PM + praM)^\ (2.17) 

then the conclusion in claim (j2.4p is equivalent to the condition that the Lie algebra gen- 
erated by the vector fields Xi, X2, and Y has full rank at each point of R^*^+^, and this 
implies that Q is hypoelliptic(See ^15j). 



Proof of Claim: Proving that the columns of {A^^vTj^o span %m amounts to showing, 
in the end, that {P<M{dl - lfai]f^Q and {P<M{dl - lfa2]f=Q span P<AfL2[0, 2tt\. They 
do so provided that all Fourier coefficients ai{m) and a2{m) of the a's are nonzero and 
that {ai{m) , ai{—m)) and {a2{m) , a2{—m)) are linearly independent, ttt, / 0. But this is 
implied by the assumption Eq. (|1.5p . This concludes the proof of the claim and hence the 
Lemma ^^. I 

Returning now to the proof of Proposition (12. 2p . we identify u{t) with the column vector 
having the four components {Q,0,ui{t),U2{t)), and we set 



u{t) = - Vf e-*^A/D,-;i X 

Xo - / ' e-'^«GM{ZM{s))ds - e^'^^^'Xi ) . (2.18) 



Putting u{t) into equation (j2.13p . we have 

XM{t) = - I e^'-'^^''GM{ZM{s))ds 
Jo 

+e*^"Xo. (2.19) 

It is then easy to see from this equation that XM^ti) = Xi. 

We consider the right side of this last equation (|2.19|) as a nonhnear operator K(Zj\,f (•)) 
which maps the Banach space C[o^4^]('Hm) (continuous functions mapping [0,ti] to T-Lm) 
into itself, with XM{t) = K{ZM){t). Let = {Z E C[o,t,](^M) : ||^(t)||oo < C} with C 
a uniform bound of the right hand side of ()2.19p . Then Q is closed and convex. Denote 
the range of the operator K by 0; then il C 0, and functions in Q are equicontinuous in t 
(Gm is bounded), so Q is compact by the Arzela-Ascoli theorem. By Schauder's fixed point 
theorem, ([23], page 151), K has a fixed point Xl,j{t) such that Xl^{t) = K(X|^)(t), which 
satisfies the conditions X|^(0) = Xq and X'lj{ti) = Xi. Here X|^ is a solution of Eq. ()2.13p 
with control n(t) given by Eq. (j2.18p . This concludes the proof of the controllability and 
hence the uniqueness of the invariant measure, Proposition (|2.2p . | 



3 Estimates on the covariance of the field in the stationary 
state 

Let ^M{x,t) = {4>M{x,t),TTM{x,t),ri^]\j{t),r2^M{t))'^ be the stationary field corresponding 
to the M-cutoff stochastic differential equation (jl.2p with invariant measure fiM- We write 

fM,n,a = ifM,n,a,4>^ fM,n,a,n^ fM,n,a,r) ^ ^ ^^ft eigenvector for Am, fM,n,a^M = Af^„ ,,, CX- 

pressing /f^^^_^ in its components, and similarly we write e^^„_^ = ifM,n,a,<p^ fM,n,a,n^ fM,n,a,r) 
for the column components of the right eigenvector corresponding to the same eigenvalue 
^Mna' "^^^ eigenvectors are normalized so that 

\JM,n,a^^M,n,a/ ('^•^) 



J dx (^f^j^^^^^^{x)e^,j^^^^^^ix) + f^j^^^^Jx)e%^^^^^Jx)^ + 



J M,n,o-,r'^M,n,(j,r 



(making e^^^ naS^ fhi naih ^^^ kernel of a projection), and, for the sake of definiteness, the 
vr-component fMnan ^^^ L^[0, 27r]-norm l/v2- For this normalization, asymptotically for 
n and M large, ej^naw ^^^ ^^^ same l/\/2 normalization. 

In this section, we obtain estimates on the covariance of <I?m(0 at equal times with 
respect to the invariant measure hm, -£'MM[^M(/M,n,CT'0 ^^(/m,™,^"^)*]' where 

^A/(/,i) = (/,^Af(t)) (3.2) 



These covariances are independent of time, and so we henceforth suppress the time t in 
their expressions. 

The variances of the rM variables are bounded by the average of the temperatures, as 
imphed by the fohowing identity. 



Lemma 3.1. We have 

E^^Jrii] ^ E^Jrl^,+rl^,] = ^{T^+T^). (3.3) 

Remark: The lemma says that the average expected energy of each of the r variables is given 
by one-half the average of the temperatures, as one might anticipate from equipartition of 
energy in equilibrium. 

Proof: Let HM{t) be the (degenerate) Liapunov function defined 

HM{t) = \\\^M{t)fH + I dxh{<t)M{x,t)) (3.4) 

^ JO 

with h{x) an antiderivative of g{x). The Ito differential of this quantity is given by 

dHM{t) = Q(ri + Ts) - rl^ dt + rMVfdoj{t). (3.5) 

By stationarity of ^m, we have that dE^j^j[H(t)] = 0, and so by the non-anticipating 
property of rMit) and thus independence of rM{t) and dw, 

= E^JdHuit)] = E^^,^[\{Ti + Ts) - rl,] dt = 0. (3.6) 



The next task is to obtain an estimate on the covariance of different modes of the 
field. In the following, \\fM,n,a,J\2 is the L2[0, 27r]-norm of fM,n,a,n^ ^nd {fM,n,a,n^^) ^^ *^6 
L^[0, 27r]-inner product of fj^.j^ ^^ with any other L^[0, 27r]-function e{x). 

Lemma 3.2. (Equal time covariance). We have, for n ^ m, or for n = m and the it 's 
differing, 

i?M„[$M(/i„,J$M(/£^,.,)*]| (3.7) 

< 7T± ±'^ T ('^[I^M(/M,„,<^)lY^^II/M,m,<x',7rll2ll5'l|oo 

+E[\^MifM,m,a')\'^]^^'^\\fM,n,a,nh\\9\\^ + {fM,n,a^'EfMX,^''^'^ 

Proof Again by stationarity and by using Eq. ()1.2p . we have that 

dE,,,mf^f,n,.Mf^lrr,.'r]=0 

= (^M,n,a + ^M,La')E^iM [^m(/m,„ „)$M(/M,m,<7')*] dt 



+{f^,n,.^Tf^iAn)dt. (3.8) 
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Ineq. ()3.7p follows immediately. ■ 

Finally, we estimate the variances of individual modes (actually their sum). We use the 
notation | e^j nar\ ~ ( ^M nar ' ^Mn ar) ' where we have made explicit the dot product 
of the two tm components for ej^^ ^^. 

Proposition 3.3. There exists a constant C > independent of the ultraviolet cutoff M 
such that 

El4,n,.,/^[l^M(/i„,JP]<C. (3.9) 

n,(T,it 

The weight \eMnar\'^ '^^ comparable to (n^ + 1) ; there exists a constant c\ with 

Clin' + 1)'-^ < \ei^n,.,f < clHn' + 1)'-^ (3.10) 

Proof. We begin by writing rj\/ = {rM,i,rM,2)'^ in an eigenfunction expansion, 

TM = Yl 4/,n,<7,r^M(/M,n,<x)' (^-H) 

n,cr,it 



(the expansion is complete!) and denoting the sum on the left side of the proposition 

|2 
\l2- 



inequality (13. 9p by ||rA,/||«2- By Lemma (13. ip above and the expansion for ru, we have that 



-(T,+T2) = E[rl,] 


= ¥m\\% 




+ y. y. 

in,a,±) ^{m,a',±') 



^M,n,(7,r ^M,m,a',r ' ^ 



(3.12) 

'^MifM*n,a)'^M{fM,m,a' 



where the double sum is over off-diagonal terms, n ^ m, or the it's different, or cr 7^ a' . 

For this equation, we use Lemma p.2p to estimate the double sum over the non-resonant 
terms with n 7^ m or the it's differing. For these terms, the denominators | A^j „ o- + ^m m a' I 
behave like |n — 7Ti| if n 7^ m and it = it', or like |n + 7Ti| if it = =p', and so in either 
case these denominators are not dangerous (see Lemma lA. II of the appendix). Thus, the 
"kernel" |A^^^ + ^Mrna'\~^ ^® £^-summable for any p > 1. We also use the fact that 
4i,n,.,r = («'eM,n,a,r)/(l + >'M,n,a) = 0{n'-^), by Eqs.(lA.2|A.6p of the appendix, which 
happens to establish Eq. (|3.10p of the Proposition as well. It follows that |ejvf no-rl ^^ ^^" 
summable for p > 1/(1 — 0). These estimates, together with Young's inequality, then show 
that the non-resonant part of the double sum in Eq. ()3.12p is bounded below by 

-Cl||rM||^2 X] l'^M,n,a,rr > -c||rM||^2 (3.13) 

\n,a,± J 

for suitable positive constants c, ci and suitably large p < 2. 
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The near-resonant terms in the double sum of Eq. (j3.12p . terms with n = m and it's the 
same but a j^ a' , are more dehcate. Fortuitously e^j ^ ^ ^ and e^ ma' r ^^^ nearly orthogonal 
for n — )• cxD, i.e., 



o±>* 



.± 



\^M,n,a;r ' ^M,n,a',r\ 



I it, * I I it, * I 

I^M,n,(T,rl \^M,n,a'r\ 



X O 



n^^'Mnn 



(3.14) 



by Lemma (jA.ip . Eq. ljA.SP of the appendix. This inequality implies that for some fixed A^ 
(independent of the ultraviolet cutoff and chosen so that the O (n^^~^ Inn)- factor is < 1/2 
for n > N), the tail series satisfies 



{{n,(7,±):n>N} 



E 



HftLnJHftlla' 



> -l^¥M\\f2. 



(3.15) 



Now the sum of the variances of the low n modes n < A^ is bounded by a constant C2, since 
the variance of an individual mode ^mUm na) ^^ certainly bounded by const/(ImAj^^ j^^)^, 
as seen from the Duhamel integral representation, Eq. (j2.2p . for $jv/ in the t — )• oo limit. 
Thus, combining Eq. (|3.12p and Ineqs. (|3.13l3.15p . we obtain a quadratic inequality for 

II^m||£2, 



:(Ti+T2)> \\rM\ 



1, 



£2 



c\\rM\\fi 



riviWfi 



C2, 



(3.16) 



which gives the bound of the proposition. 



In the inequality (j3.9p of Proposition (13. 3p we can actually replace function fj^ no-™ ^^^ 

expectation E |*&(/jvf „ o-)l^ with the free eigenfunctions fMna^^)^ ^^^^ ^^' ^^^^ eigenfunc- 
tions of Bj\,f , i-e. the matrix operator obtained from Am by setting the coupling functions a 
to zero. For n > 0, we take /a/^„_^(0) = -^{±i{n^ + l)-^/'^e'''"-'' ,€'""■'' ,0), with cr = 1 or -1 
(except when n = 0, where there is no a-dependence) , and /m,-i,(t(0) = (0,0,1,0) or 
(0,0,0,1). 

Corollary 3.4. There exists a finite constant C > independent of the cutoff M such that 



n,a,± 



< c. 



Remark: The corollary implies that 



Y,{n'' + l)'EUM{n)\^]<C 



(3.17) 



(3.18) 



uniformly in M, with (j)M{n) = j e~'^'^^(j)M{x) the usual Fourier coefficient of the field com- 
ponent (l)M of <l>Af. If > 0, then the field 4)m is in L?' a.s., with the variance of ||i;^Af lb 
uniformly bounded in M; obviously, each mode 4>Min) has variance uniformly bounded in 
n. However, we believe that this bound can be improved, that in fact var((/)Af (n)) = 0{n^'^) 
and that var(7rA//(ra)) = 0(1) uniformly in M as in the linear and non- linear equilibrium 
cases and the linear 5 = non-equilibrium case. This improvement will be the subject of a 
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subsequent investigation. 
Proof: We expand 



^MU^,n,M)= E (/A/,n,.(0)'4r.,.')^M(/f;,^,.0- (3-19) 

m,a' ,±' 



The coefficients of this expansion (with inner products as in Eq. (j3.ip ) satisfy 



0{m '^), n = — 1,?TT. — )■ oo 
0{'nP~^), m = —l,n — t- oo 

0(j^), m/n, m,n^oo 

0{1), otherwise, 



(3.20) 



ah uniform in the cutoff M (see Eq. ()A.2ip where these estimates are shown). Substituting 
the expansion Eq. (j3.19p into EH^mUm n (t(0))P]5 ^^ obtain the double sum of terms 

Off-diagonal terms m ^ m' must be estimated using the non-resonant Ineq. ()3.7p of Lemma 
(|3.2p and the above inequalities (|3.2Up . For m,m' ^- oo these terms are 

^ 7 V TVf ^^ 3.22 

\(n — m)[n — m')(m — m) J 

It remains to sum the on- and off-diagonal terms ()3.2ip over m,m', and over n with 
the prefactor n^^~^. One uses the fact that fn^~^E[\^M{fej^ ^^il"^]^'"^ is £^-summable in 
m, the content of Proposition (|3.3p . The near-resonant m = m' terms are estimated using 
the inequalities (|3.20p . By extensive use of Young's inequality and Holder's inequality, the 
triple sum is shown to be finite. ■ 



4 On the tightness of the ultraviolet cutoff stationary mea- 
sures 

Let {fJ-Ai} denote the unique stationary measures for the ultraviolet cutoff systems of the 
previous sections, M labeling the cutoff. Let ^uit) be the canonical stationary process 
associated with fiM so that in particular hm is the law for $Af(i) for any time t. It will 
be convenient to regard the field ^M^t) as taking values in a space of distributions dual 
to the Schwartz space S = C^^[0, 27r] x C^JO, 27r] x M x M, but where the Fourier modes 

{(puin), f^M{n)} of $Af are all zero for \n\ > M. 

Corollary ()3.4p provides a bound on the variance of the Fourier mode of the field, 
^;ar^J^^(<I>^f (/j^^(O)) ~ n^~^^ uniform in M at a fixed time, say t = 0. As before, we 
suppress the explicit time dependence. Let Tig = Hg (B Hg-i ® C ® C. Then generally, we 
note the following: 
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Proposition 4.1. Assume that the fields ^m <ire of mean zero with respect to /xj\,/ for all 
ultraviolet cutoffs M , and that the field modes are of variance 



var 



^,,n^M{ftAm^E,,A\^M{ftAm^]<C{l + n'Y (4.1) 

for some constants C and power p independent of M . Then the measures {/xjv/} are tight 
in the weak-* sense that there exist a subsequence {fJ-Mj}, j = Ij 2, ..., a limiting measure fi, 
and a limiting field $ such that 

lim E [F{^M, (/l), $M, (/2), ..., ^M, Uk))] 

= E^[F{^fi),^{f2)....,Hfk))\ (4.2) 

for all hounded continuous functions F{xi,X2, • • • , x^) on R , and with /i, /2, ...fk G S for 
all k. The limiting measure has support in the space of distributions S' and is a -additive 
on the Borel a-algebra generated by cylinder sets of the form 

{^ : ($A/(/i), ^Af (/2), ..., ^M{fk)) G B} (4.3) 

with base /i, /2, -.fk G S and B a Borel set in Mr. 
We have that for s < —p — 1/2 

E^mWns] < Ci (4.4) 

for some finite constant Ci; the same bound holds for the hm 's- In particular, ^ is in %s 
a.s. 

Remark: Again, if the Fourier coefficients of the a's behave as a power, |a(n)| ~ n , then 
the variances of the Fourier modes indeed satisfy Ineq. (|4.ip . and so we obtain a hmiting 
measure for which ^ € T-Lg^ s < 9 — 3/2. However, we do not claim that this measure is an 
invariant measure for the nonhnear stochastic wave equation; this remains an open problem. 

Proof: The marginals of {/^Af } restricted to functions F just depending on a fixed and finite 
number of Fourier modes {^M{fn,ai^))}{n: \n\<no} ^-re tight. This is the case since a closed 
ball in a finite-dimensional Euclidean space is compact, and given e > 0, 

P,A'^:\\P{ln\<n.}n2>R} < ^E^,A E I^A/(/„^.(0))P] 

|n|<no,CT,it 

^ :^ E i^ + n'r<e (4.5) 

|n|<7io,CT,it 

for another constant c for R sufficiently large, by our bounds on the Fourier coefficient 
variances (see [28j). Here, P||„|<„^}<1> is projection of <I> onto its Fourier modes with \n\ < Uq- 
By this tightness, one can construct a subsequence {^m „ } with convergent marginals based 
on the random variables {^M,{ftA^))}{n: |n|<no}- 

One then passes to a sub-subsequence to get convergence for marginals based on a 
larger collection {$Af(/jta-(0))}' I"-! — "-i ^^ "-i > ^o and diagonalization, one obtains a 
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subsequence {fJ-M } which, integrated against any continuous function F of a finite number 
of Fourier modes, converges, i.e., linij^oo E^j^j_{F) exists. 
Again by Chebyshev, given e 



^mm{^ ■■ min-r'-^ > R} 



< ^E 

— j^2 t^M 



rtjCjit 



< 



C ^ (l + n2)P 



n,(T lb 



for R sufficiently large, uniformly in M. This implies the last inequality of the lemma, 
Ineq.(j43D. 

It is then easy to see that the domain of definition of \x extends uniquely to functions of 
the form F($(/i), $(/2), ..., $(/a:)), with fij2,-,fk G 5, F{xi,X2, ...,Xk) bounded contin- 
uous. Ineq. (j4.6p reduces the problem of showing the convergence of {E^j^j [F]} to showing 
that of {E^^j [F; -B/j]}; where Br = {<!> : ||<5||-Hs < R}- We write f<N as the projection of / 
onto the subspace of Fourier modes with \n\ < N and />Ar for the tail of its series. Now for 
^ e Br and for any 6 > 0, \{fj, ^) - {fj,<N, ^)| < R\\fj,>N\\n^s < (^ for iV sufficiently large. 
The uniform continuity of F on the bounded set {x : \xj\ < R} then gives the convergence. 

That fj, is countably additive on the Borel sets generated by the cylinder sets is also a 
consequence of Ineq. (l4.6p ([13]). This concludes the proof of the proposition. | 



A Perturbation theory for the eigenf unctions of the operator 
A 

This appendix summarizes properties of the left eigenfunctions {/„ o-} ^^'^ right eigenfunc- 
tions {e^g.} and their corresponding eigenvalues {A^^.} for the matrix operator A. All 
estimates are uniform with respect to the ultraviolet cutoff M, and so the index M is sup- 
pressed. We assume throughout this appendix that the coupling functions ai and 02 satisfy 
the assumptions (II. ip . See ^16j for the perturbation theory methods utilized. 

Fixing the mode number n > and the it's and then suppressing these and other indices 
except as needed, we write /^^ = {fa,4>, fa,Tr, fa,ri, fa,r2) foi' the components of f^a: ^^'^ ^^ 
write e^^ = {6^,4,, 6^,77, ea,ri, fa,r2)'^ for the column components of the right eigenfunction. 
The corresponding eigenvalue A„ ^ is near iti vr? + l ; for n > and for a given choice of 
lb, the eigenvalue is nearly doubly degenerate, whence the index o" = 1 or — 1. 

From the eigenvalue equation {f^a > -^ ')'H — ^naifna > ')'H^ o^^ finds the relations 



u^ = ^f.Adl-i), u= f:';'"\ (A.i) 

Act i + Ao- 



and similarly, from Ae^^ = A^^-C^ 



n,cr ''n,cT n,a^ 

1 (a,e<j,7r) 



Ao- i + Ao- 



15 



Note that fo^r and eo-,r each have two components, e.g., f^^r = (/o-.n, /o-,r2)) corresponding 
to the two components of a in the above relations. One can then use these relations to write 
the eigenvalue equation as a non-linear eigenvalue equation just involving the vr-components, 

XlU. = faAdl - 1) - T^{fa,.,a){a\ (A.3) 

i + Act 

or 

A^e^,^ = {dl - l)ea,n - " \a){a, e^,^). (A.4) 

i + Act 

In these equations, |a)(a| = X^i=i |aj)(aj|- For large n, Ao-/(l + Ao-) = 1 ibi/n + ©(n^^), so 
that these latter two equations are nearly self-adjoint eigenvalue equations for 

A^e^,7r = {dl - l)e^,,r - lia{a, e^,^), (A.5) 

with Ao- determined implicitly by substituting in /? = Ao-/(l + Ao-)- Evidently, e^^^^ = 
— /3(A^ — dl + l)~-^a(a, eo-,7r), and from this representation, one sees that the m*^ Fourier 
coefficient eo-,7r(^Ti) of ecr,7r, is 0(-tz — j), for rn?' ^ m?. This implies in particular that 

(Q,e,,,) = 0(n^). (A.6) 

In the following, let Pq = -Po(^) be the projection onto the eigenspace spanned by e*"^ 
ande"™^ inL2[0,27r]. 

Lemma A.l. We have that for n large, the eigenvalues {A^^} of A are given by 

Xt = ±('v'^?^+^ii|^) (A.7) 

+ iO{n*^~^ Inn) + 0{n'^^~^ Inn) 

where the Hn^a CLfe the two eigenvalues of the operator —PQ{n)a){aPQ{n), and where the 
error terms are, respectively, imaginary and real. 

Let e^o-r ^6 ^he r-components of the right eigenvector e^^ of A. Then for a / a', we 
have 

± * i 

■ = 0(n mn), (A. 8) 



where the dot indicates the dot product of the components, \ ■ \ being the Euclidean length. 

Sketch of proof. For n > 0, let P = P(n) be the projection onto the subspace spanned by 
the two eigenvectors corresponding to the two eigenvalues A^ near — (n^ + 1) for the operator 
d^ — 1 — /3\a){a\. {P is close to Pq.) Then the shift in the eigenvalues is determined from 
the 2x2 matrix equation 

(A^ + (n^ + l))PoPPo^^ = -(3Poa){aPPoC^ (A.9) 

for ^o- in the span of Pq. The projection P itself can be estimated from its representation as 
a contour integral of the resolvent [d"^ — 1 — /3|Q;)(a;| — z)~^ expanded in a Neumann series. 
One finds that 

Xl = -(n2 + 1) + /3;u, + (3^0{n^'-' Inn), (A.IO) 
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where fia is one of the two eigenvalues of the rank 2 matrix —Poa) {aPo and is of the order 
of |a(n)p = sup |dj(n)p = 0{'n?^). The Neumann series for the resolvent above is in powers 
of (a, (5^ — 1 — z)~^a) with z traversing a circle {z : [z + n^ + 1| = n/2}. These powers are 
estimated using 



a[m]\^ v-^ |a;(n)P 



I (a, (9^ — 1 — z) ^a 

{m:m<n/2\ {m,:n/2<m<3n/2,m7^ra} 

|Q;(n)p Y^ |Q(m)p 



,9 91 



{m:r?i>3n/2} 



n , ■^^ , , m? 



< c|a(n)plnn ^ 

~ n 

(The second sum in the second line accounts for the In n/n factor, by an integral test.) One 
uses this estimate to obtain the correction to the eigenvalue shift in Eq. ijA.lOp . See [30] for 
additional details. We then set /3 = Ao-/(l + Ao-) into Eq. (|A.10p to determine Ao- implicitly; 
this gives the first assertion of the lemma, Eq. (|A.7p . 

We also need an approximate orthogonality relation for ea,r and e^'^n cr ^ cr' ■ In 
Eq. ()A.4p . we decompose eo-^jr as 

where ^^ is an L^-normalized eigenfunction of PQa){aPo, and hence in the subspace spanned 
by Pq. Then 

Po{\l + n2 + 1 - ^a){a\){C + P^e^,. + Qoe.,^) = 0, (A.13) 

where Qo = 1 ~ ^o- We have that 

\\PQ{\l+n^ + l-pa){a\)C\\ =0{n^^'^\nn) (A.14) 

by the eigenvalue shift estimate above, Ineq. (|A.10p . Also, we have that 

Qoea,^ = -PQo{\l-dl + l)-^a){a,e^^^), (A.15) 

so 

|(aQoe.,.>| = 0(n=^'-'lnn), (A.16) 

and thus 

||Po(A^ + r? + l- /3a)(a|)Qoe.,OII = 0{n^'-^ Inn). (A.17) 

From the identity Eq. ijA.lSp . this last equation, and Eq. (|A.14p above, it follows that 

{n„ - ^JL„, + 0{n^'^-^\nn)\ \\Poe„^A = 0(n^^"^ Inn), (A.18) 

where a' refers to the complementary value of a. The assumption on the coupling functions 
a, Eq. ljl.ip . assures that \^fj — fi^i\ > c|d(n)p, so that 

\\Poe^,n\\=0{n'^'^-^lnn). (A.19) 
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The last equation of the lemma, Eq. (jA.SP , follows from the equation for eo-,r ) Eq. ()A.2p ; 
the decomposition for eo-.Tr in Eq. ()A.12P ; and the orthogonality of ^^ and ^^, under Poa) {oPq. 
We have that 

= (C°,")(",Qoea',7r) + {QoecT,7T,a){a,C',7T) + (^oc^.tt, ")(", Qoe<x',7r> 

+ (Qoea,7r,a)(a,-Poecr',7r) + (Qoe^.Tr, «)(«, Qoea'.Tr) 
= O(n^'^-Mnn) (A.20) 

by Eqs. (|A.16|A.19p . Combining this equation with Eq. (jA.2p . we obtain Eq. (|A.8p of the 
lemma ■ 

Eq. (j3.20p of the text provides estimates on inner products (/^o-(O), e^ o-')'^- '^° illustrate 
how these estimates are obtained, consider a case m ^ n, m,n large. By familiar resolvent 
identities, and with B the unperturbed matrix operator with the a's turned off and with 
the left eigenfunction fnai^)^ 

(&(0).<,.>« = {/*,(0),Fo(n)(F(m) - -Po(m))4' .)« 



'' Ym. 



27ri .,^, 



e± /'n^ ^,\o±' 



1 f ^ {f^^^JO),a)e^^^, 
az 



,ty ,< 



2vr^7,„ (A±(0)-z)(A±;,-z) 

a{n)a{m) -q flfnf^] (a 21) 



(A^(o)-A^;,„,.,)(A^;,™,.. + i) V^-"^ 

Here, 7m is the circle {z : \z — A,^^, | = 1/2} and, in the second line, C = A — B is the 
4x4 perturbation matrix of just the a entries, all other entries being zero. In the last line 
we have used Eq. (IA.4[) for e^^, ^, and the fact that (a, e^ ^, ^) ~ a{n) = 0{n^), Eq. (lA.6p . 
The other relations of Eq. (j3.20p are analyzed similarly. 
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